Cronin effect vs. geometrical shadowing in d+Au collisions at RHIC 
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Multiple initial state parton interactions in p{d) + Au collisions are calculated in a Glauber- Eikonal 
formalism. The convolution of perturbative QCD parton-nucleon cross sections predicts naturally 
the competing pattern of low-pr suppression due geometrical shadowing, and a moderate-pr Cronin 
enhancement of hadron spectra. The formal equivalence to recent classical Yang-Mills calculations 
is demonstrated, but our approach is shown to be more general in the large x > 0.01 domain because 
it automatically incorporates the finite kinematic constraints of both quark and gluon processes in 
the fragmentation regions, and accounts for the observed spectra in elementary pp ttX processes 
in the RHIC energy range, -^s ~ 20 — 200 GeV. The Glauber-Eikonal formalism can be used as 
a baseline to extract the magnitude of dynamical shadowing effects from the experimental data at 
differente centralities and pseudo-rapidities. 

PACS numbers: 12.38.Mh; 24.85.-|-p; 25.75.-q 



It is well known that in proton (p), or deuteron (d), 
reactions involving heavy nuclei {A ~ 200) at y/s < 40 
AGeV, the moderate transverse momentum {pT ~ 2 — 6 
GeV) spectra are enhanced relative to linear extrapola- 
tion from p + p reactions. This Cronin effect Q, |^ 
is generally attributed to multiple scattering of pro- 
jectile partons propagating through the target nucleus. 
The data can be well accounted for phenomenologically 
by adding a random Gaussian transverse kick 5k^ cx 
u^L/A to the projectile partons prior to hadronization 
^ la, 0, 01 • Here A is the parton mean free path in the 
nucleus, L oc A^^^ is the average path length, and /i is 
a typical screening mass in ground state nuclei. These 
models naturally predict a slowly decreasing Cronin ef- 
fect with increasing energy which has only recently been 
possible to test at ^/s = 200 AGeV at the Relativistic 
Heavy Ion Collider (RHIC). 

Interest in the Cronin effect has been revived, due to 
the developement of a new formulation of the physics 
based on the concept of gluo n saturation and classical 
Yang-Mills field models In addition, a radical pos- 

sibility was proposed in Ref. [Sj , that nonlinear gluon sat- 
uration may in fact strongly suppress moderate-pT spec- 
tra at RHIC. If such a deep gluon shadowing in this kine- 
matic range were true, then an anti-Cronin suppression 
should have been observed in the large x — 2pT/y/s « 
0.01 — 0.1 and moderate 1 < Q = px < W GeV scale 
ran ge a ccessible at RHIC. Four experiments at RHIC 

mill El El found independently that the nuclear 

modification factor, RdAiPx) = '^.dadAl Adapp, was con- 
sistent with a positive Cronin enhancement of hadrons 
with 2 < pt < 5 GeV. The magnitude of the enhance- 
ment is somewhat smaller than predicted for the tt", as 



reviewed in jl5| , but no evidence of strong shadowing was 
reported. 

After the release of the RHIC data, the satura- 
tion model predictions of suppression were revised in 
Ref. El El 

Il8l |. In Ref. Ell Cronin enhancement with 
no high-pT suppression was shown to be a generic fea- 
ture of saturation models. In Ref. jlJi] another version of 
saturation dynamics with Cronin enhancement coupled 
with the high-pT suppression of Ref. was discussed. In 
Ref. Ell I the Cronin enhancement at y = was predicted 
to be progressively negated by non linear QCD evolution 
at smaller nuclear x, and therefore a gluon shadowing 
suppression is predicted at higher rapidities. 

Different approaches to the calculation of the Cronin 
effect can be formulated in infinite momentum and tar- 
get frames. In the traditional Glauber-Eikonal (GE) ap- 
proach [l^[2olEl| . sequential multiple partonic collisions 
in the target frame are computed. This leads to trans- 
verse diffusion and unitarity is naturally preserved. The 
\ow-pT spectra are suppressed by unitarity to compen- 
sate for the moderate px Cronin enhancement. This is 
what we call "geometrical shadowing" , as it is driven by 
the geometry of the collision. No high-py shadowing is 
predicted in this approach. 

In applications, the GE series has been directly evalu- 
ated thusfar only up to the three-scattering term and for 
\/s < 40 GeV p, 'l!/|. Numerically more convenient ap- 
proximated GE models 4, ^ H Q have been proposed. 
They modify the pQCD rates through the inclusion of 
a nuclear broadened intrinsic fcx, instead of evaluating 
the full GE series. Phenomenologically, it is well known 
[23. that intrinsic fc^ ~ 1 GeV must be introduced to 
correct coUinear factorized pQCD predictions to account 
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for p + p data at moderate pT < 5 GeV. The approxi- 
mated GE models simply extend that idea by adding a 
random kick 6k^ to the intrinsic k^. One drawback of 
such approaches is that a non-trivial px or collision num- 
ber dependence of the effective nuclear transport coeffi- 
cient AtVA ~ 0.05 GeVVfm Q must also be introduced 
to account for the actual Cronin data. While a logarith- 
mic pt dependence of Sk^ipx) is expected for partons 
undergoing multiple Yukawa screened interactions [2l|. 
the functional form of that px dependence is usually ad- 
justed to fit the Cronin data at one energy. A further 
drawback of such approximated GE models is that the 
unitarity constraints built into GE are ignored and hence 
the unitarity shadow and Cronin are treated as two seper- 
ate phenomena. 

The more recent approaches [13, IH, H^l to the 
Cronin effect in the infinite momentum frame are based 
on the Mclerran-Venugopalan (MV) model of the nuclear 
wave funtion in classical Yang-Mills theory [i^ . The gen- 
eral equivalence of GE and MV formulations for trans- 
verse diffusion was discussed in [sO, |^ in the context of 
gluon dominated small x 1 kinematics. In these ap- 
proaches, the nucleus is approximated by a Weiszacker- 
Williams gluon field with non-linearities approximated 
semi- analytically or computed numerically |27l l29l|. The 
non-linear gluon interactions lead to transverse diffusion 
and hence Cronin enhancement of nuclear partons prior 
to the scattering. The essential scale in this approach is 
a gluon saturation scale Qs = Qs{y,y/s,A), with y the 
rapidity of the produced gluon. One of the advantages 
of the MV approach is that unitarity is at least approxi- 
mately enforced through the conservation of the number 
of virtual gluons in the transverse diffusion. Therefore 
these models predict a definite anti-Cronin suppression 
below some scale oc Qg. On the other hand, a disadvan- 
tage in present formulations of the MV model is that they 
ignore finite-energy kinematics of valence- and sea quark- 
induced processes and the non-asymptotic large x > 0.01 
features of gluon structure, where the classical approxi- 
mation is unrealistic. A major disadvantage of MV mod- 
els is that they cannot account for the elementary p -\- p 
transverse spectrum, that forms the denominator of the 
RpA nuclear modification factor. Neither can the models 
reproduce the absolute normalization of the spectra in pA 
collisions without extra phenomenological assumptions. 

In this letter, we compute directly the GE series via 
numerical convolution of elementary parton-nucleon pro- 
cesses. An advantage over approximated GE models 
0, 0, 0, 0| is that our approach automatically conserves 
unitarity in the geometric optics sense of GE theory. 
In addition, we do not introduce extra phenomenolog- 



ical pT-dependent nuclear broadening of the intrinsic 
kx, since the GE series predicts the functional form of 
the Cronin enhancement based on the calculated non- 
asymptotic pQCD parton-nucleon cross sections. An ad- 
vantage of our GE approach over MV model applications 
is that we autmatically include the finite kinematic large- 
X feautures of both quark and gluon processes. Perhaps 
the most important advantage over the MV approaches 
is that our formulation is directly constrained to repro- 
duce the absolute normalized spectra in p + p collisions. 
Therefore the GE approach presented below calculates 
consistently both p + p and p + A collisions at the finite 
energies accessible at RHIC. 

Beside the geometrical quark and gluon shadowing, 
which is automatically included in GE models, at low 
enough x one expects genuine dynamical shadowing due 
to non-linear gluon interactions, as described in satura- 
tion models. Both kind of shadowing are present in the 
data, but it is not possible a priori to tell in which propro- 
tion. The Glauber-Eikonal formalism can then be used 
as a baseline to extract the magnitude of dynamical shad- 
owing effects from the experimental data. 



PARTON-NUCLEUS COLLISIONS IN THE 
GLAUBER-EIKONAL MODEL 

The GE expression for a parton nucleus scattering |2£ 
El is: 



da 



iA 



cPpxdycPb 



CO „ 
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n=l 

da'^ 



da'^ 



X e 



-^■"(Po)Ta(6) 



S{Y.h-pT), (1) 



where 7^(6) is the target nucleus thickness function at 
impact parameter b. The differential and integrated 
parton-nucleon cross sections, da"^^ /d'^k and a'^^ {po) — 
J da"^^ , are computed in pQCD as discussed below, see 
Eq. lfSa|) - (|5b|l . The integrated parton-nucleon cross sec- 
tion depends on an infrared scale po, which is determined 
by fitting p + p data. The exponential factor in Eq. 
represents the probability that the parton suffered no 
semihard scatterings after the n-th one. In such a way, 
unitarity is explicitly implemented at the nuclear level, 
as discussed in Ref. The sum over n starts from 

n = 1 because we are interested in partons which are 
put on-shell by the interaction and later on hadronize. 
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The sum over n may be performed in Fourier space. The 
result reads: 



da 



iA 



(PpTdycPb J An^ 



e-^'="-^^-5'^(rT;po) , (2a) 



where, suppressing the dependence of the quantities in 
the r.h.s. on y, 

S'^{rT;po) ^ e--'"('--Po)T.(6) _ g-.-(po)T.W (2b) 



and 



iN 



irT;Po) 



-ik-TT 



da 



iN 



d^kdyi 



(2c) 



It is possible to show [20j that, in the high-pT- limit, 
Eq. O reduces to the usual single scattering approxi- 
mation for parton-nucleus scattering: 



da' 



iX 



dp^dyd^b 



Pt- 



TA{h) 



da 



iN^iX 



dpi^dy 



(3) 



As pt^O unitarity corrections switch on, suppressing the 
integrated parton yield [s^ , and inducing a random walk 
of the parton in px space, thus redistributing the partons 
to higher pr compared to the single scattering approxi- 
mation j2i|. This is how the multiple scattering mecha- 
nism of Eq. induces "geometrical" shadowing at low 
Pt, and Cronin enhancement of the transverse spectrum 
at moderate pt, respectively. Note also that the accu- 
mulation of kr kicks in the multiple scattering process 
is computed in the model, not input as a Gaussian fold- 
ing as in approximated GE models. The full expression 
for the PT spectrum, Eqs. and if^ - lf^ . interpolates 
naturally between the geometrical shadowed low-p^ and 
the Cronin enhanced moderate-pT regions. 

Note that a^^{rT) (x as tt ^ and a"^^ {tt) 
(T*^ as tt oo. This suggests the interpretation of 
o'^^ [tt) as a dipole-nucleon "hard" cross section. This 
dipole is of mathematical origin, and comes from the 
square of the scattering amplitude written in the Fourier 
variable r^, which represents the transverse size of the 
dipole. Then, we can interpret 5*^, Eq. (|2bp . as the 
dipole-nucleus "hard" cross section, which clearly incor- 
porates Glauber-Gribov multiple scatterings of the colour 
dipole. No other nuclear effects on PDF's are included 
beside multiple scatterings. 

The interpretation of the interaction in terms of multi- 
ple scatterings of a dipole allows to relate this approach 
to other multiscattering formalisms such as Ref. ISTl ISS l 
and the saturation computations of Ref. 
A first step in building such a dictionary was taken in 



Ref. [321, where it was shown that the dipole cross sec- 
tions of Eqs. (|2b|) - (|^ are equivalent, in a suitable kine- 
matic region, to the dipole cross section considered in 
the MV model of Ref. Hi H^]. The main difference is 
in the input parton-nucleon cross section in pp-space. 
In our case, as we will discuss in the next sections, it 
is computed in pQCD, including full kinematics and in- 
teractions of the incoming parton with both quarks and 
gluons. Moreover, its energy and rapidity dependence 
are controlled by the DGLAP evolution of the parton 
distribution functions of the target nucleon. In the MV 
model, it is approximated using asymptotic kinematics 
for gluon targets only. Both models consider in 5'"* only 
inelastic dipole-nucleus scatterings. They neglect diffrec- 
tive dipole-nucleus interactions, which however modify 
the pT-spectrum only at the lowest tarnsverse momenta 

Beside the geometrical quark and gluon shadowing, 
which is automatically included in GE models, at low 
enough x one expects genuine dynamical shadowing due 
to non-linear gluon interactions as described in the satu- 
ration models. However, it is difficult to disentangle these 
two sources of shadowing and suppression. The distinc- 
tion between the two is however of fundamental interest 
as has already been emphasized in the context of e -I- p 
DIS HERA by A. Caldwell H. Most theoretical interest 
is not in the ubiquitous geometrical shadowing and uni- 
tarity corrections, but in the onset of genuine nonlinear 
QCD physics ^1|. Moreover, saturation models cannot 
predict as yet the upper bound on x below which non lin- 
ear effects set in. In order to help recognize possible novel 
nonlinear regimes it is essential to be able to calculate the 
baseline spectra isolating the unitarity and geometrical 
shadowing alone. The parton level GE model discussed 
below provides such a baseline. 



INCLUSIVE MINUET AND HADRON 
PRODUCTION IN pp COLLISIONS 



Let's consider a pp' collision, where p and p' stand for 
a proton (p), a deuteron (d), or a nucleon {N). In leading 
order pQCD, the inclusive cross section for production of 
a parton of flavour i — g,q,q {q = u, <Ls, . . . ) with trans- 
verse momentum pT and rapidity y [221 may be written 
as a sum of contributions to the cross section coming 
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from projectile (p) partons and from target (p') partons: 



da 



pp — nX 



dp\dy 




da'P 



'-"^ dy.d^PT 
da'P 



yi=v 



dyid^PT 



Here we considered only elastic parton-parton subpro- 
cesses, which contribute to more than 98% of the cross 
section at midrapidity [s^l- In Eq. Q), 



(4) 



{-^fi/p)yi,PT 

da'P' 
d'^prdyi 



1 f d& 
1 4, j ^y2Xifi/pixi,Ql)--^ {s,i,u) 



^■2fj/p' ix2, Ql) 



da'P 
d^prdy.^ 



da - 

dy2-^ {S,t,u) Xljjlpi {X2,QI) 



(5a) 
(5b) 



are interpreted, respectively, as the average flux of in- 
coming partons of flavour i from the hadron p, and the 
cross section for the parton-hadron scattering. The ra- 
pidities of the i and j partons in the final state are la- 
belled by yi and y2- Infrared regularization is performed 
by adding a small mass to the gluon propagator and 
defining niT = \/p't + Po ■ fractional momenta of 

the colliding partons i and j are 0:1^2 — '^(6*^' -|-e*^^), 
i.e. the incoming partons are coUinear with the beams. 
The integration region for j/2 is — log{yG /rriT — e~^') < 
Z/2 < log('\/s/"T^T — e^'). The summation runs over parton 
flavours j = g, q, q. The partonic Mandelstam variables 
are 

i = -m|(l + e"^'+^^) 

S = —t — U = X1X2S . 

For the parton distributions we use the CTEQ5 
parametrization at leading order The choice of the 
factorization scale Qp is discussed later. The cross sec- 
tions da^^ /dt of the ij—ni elastic partonic subprocesses 
can be found, e.g., in ^3. They are proportional to 
as(M^), computed as in |39|, at a scale /i = Qp The fac- 
tor K in Eq. is introduced in order to account for 
next-to-leading order (NLO) corrections 0|i is in 
general \/s and scale dependent 0,113. 

Inclusive hadron production through independent frag- 
mentation of the parton i into a hadron /i, is computed 
as a convolution of the partonic cross section (0J with a 



fragmentation function Di^h{z, Ql)- 

dcjPP'^'^^ daPP'^"^^ 
dq^dyh dpj.dy^ 

where qx is the transverse momentum of the hadron h, 
yh its rapidity, and z the light-cone fractional momentum 
of the hadron and of its parent parton i. For details, 
see Eqs. (8)-(ll) of Ref. 39]. In this letter, we use LO 
Kniehl-Kramer-Potter fragmentation functions 01 and 
set the fragmentation scale Qh — Qp- 

In the computation of the pp' cross section we 
have two free parameters, po and K, and a somewhat ar- 
bitrary choice of the factorization, renormalization and 
fragmentation scales. Our strategy is to compare two 
choices for those scales, namely Qp = Qh = mr/^ and 
Qp = Qh = TnT, and then fit poi to hadron production 
data in pp collisions at the energy of interest. We ana- 
lyze here tt^ production at -/s = 27.4 GeV 0, and 7r° 
production at ^/s = 200 GeV 0|. For the i^- factor we 
perform a fit to the high-pr tail of the data, following 
the procedure described in Ref. [s^ • The fit of pq is per- 
formed by requiring that the computed spectrum does 
not exceed the experimental data at low pt ^ 1 GeV. 
At = 200 GeV, this fit is difficult because data exist 
for qT > 1.2 GeV only, so we used also data on charged 
hadron production |33l |. The resulting data/theory ratio 
is plotted with thin lines in the bottom panel of Fig. ^ 
and the extracted parameters are listed in Tabled Note 
that the if-factor is strongly correlated to the choice of 
scale, while po is more stable. Both of them depend on 
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(4> 


Qp — Qh 


y/s = 27.4 GeV 


= 200 GeV 


0.52 GeV^ 


m.T/2 


po = 0.70 ±0.1 GeV 
K = 1.07 ±0.02 


PO = 1.0 ± 0.1 GeV 
K = 0.99 ± 0.03 


mj- 


po =0.85 ±0.1 GeV 
K = 4.01 ± 0.08 


po = 1.2 ±0.1 GeV 
i^ = 2.04 ±0.12 


GeV^ 


mT/2 


PO = — 
K = 3.96 ± 0.11 


PO = — 
K = 1.04 ± 0.06 


rriT 


Po= — 
K = 13.4 ±0.4 


PO = — 
ii: = 2.04 ±0.12 



TABLE I: Fitted values of po and the /('-factor for tt and 
7r° production in pp collisions at ^/s = 27.4 GeV and ^/s — 
200 GeV, respectively. We quoted the fit uncertainties only. 
The systematic uncertainty in the absolute normalization of 
experimental data (20% and 9.6%, respectively), which affect 
the determination of the K-factor are not included. The fit 
of the if-factor in the case of zero intrinsic momentum was 
made with the po determined using the (fcy) = 0.52 GeV'^. 



As in Ref . |39j , we obtain a satisfactory description of 
data for (/t ^ 5 GeV over a broad range ^/s, but the 
curvature of the hadron spectrum is overpredicted in the 
qt = 1 — 5 GeV range. As it is well known jij lisj . 
this can be corrected by considering an intrisic trans- 
verse momentum kx for the coUiding partons '42] . There 
exists many ways of implementing it phenomenologically, 
and we choose for simphcity a kx smearing of the cross 
section to approximate this effect. We introduce then 
unintegrated parton distributions 



fi{x,kT,Ql) 



Tr{k, 



2 \ fi ("^i Qp) 



where the width {k'^) of the Gaussian enters as a phe- 
nomenological parameter, and convolute over cPkiT and 



d k2T in Eqs. l(5ajl - l(5b(l . 

We found that a fixed {k'^) = 0.52 GeV leads to a dra- 
matic improvement in the computation of the trasverse 
spectra, which now agrees with data at the ±40% level. 
The quality of our pQCD computation including intrinsic 
kx is shown in Fig. ^ and the extracted K is reported 
in Table HJ Without intrinsic kr it is not possible to fit 
the value of po, due to the steepness of the data/theory 
ratio. The fit of the iiT-factor in this case was made with 
the Po determined using the intrinsic /ct- The optimal 
choice of scale is found to be Qp = Qh = mT/2 at both 
energies, as the value of the iC-factor is the closest to 1. 

Let us comment briefiy on the physical meaning of 
the infrared regulator po- The divergence of the pQCD 
cross section for minijet production, indicates clearly a 
break-down of unitarity at low pt- A phenomenological 
way of restoring it, is to tame the divergence of da^^ /di 
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FIG. 1: Top panel: pion transverse momentum spectrum at 
= 27.4 GeV and ^ = 200 GeV. Solid lines are LO pQCD 
computations according to Eq. J^, with (fcy) = 0.52 GeV^ and 
Qp = Qh = mT/2. The regulator po and the /C-factor are given 
in Table Bottom panels: the data to theory ratio for different 
choices of parameters. Solid lines are for Qp = = mT /'i, dashed 
lines Qp = Qh = my. The pair of thin lines is computed with no fcy 
smearmg, (fc2^) = GeV^, and the pai r of thick lines is computed 
with ) = 0.52 GeV^. In the tt'' case 1351 , the dashed area shows 
the relative statistical and point-to-point systematic error added 
in quadrature, and does not include the systematic uncertainty of 
9.6%]On the absolute normalization of the spectrum. In the 
case 12, the shaded area includes statistical error only, without a 
systematic uncertainty of 20% on the absolute normalization of the 
spectrum. 



by adding a small mass regulator, po, to the exchanged 
transverse momentm pT- Seen in this fight, po repre- 
sents the scale at which higher order parton processes 
enter into the game, beside the single scatterings consid- 
ered in Eqs. (|5a|l - (|5b|l . As the center of mass energy is 
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increased, the partons probe the nucleon at smaller x, 
so that the density of target parton increases, and one 
should expect unitarity effects to arise at larger scales. 
The slight increase of the fitted po = PoiV^) with energy 
is indeed a consistent check of this picture. For the same 
reason, and since quarks interact more weakly than glu- 
ons, one might expect also the unitarity corrections for 
quark-nucleon scattering to arise at Po|quarks ^ Polgiuons- 
However, to check this relationship, we would need at 
least data on production in a pr range of 1 — 6 GeV. 
Here we simply set po|quarks = Po|giuons = Po- 



FROM pp TO pA COLLISIONS 

Having fixed the free parameters in pp collisions, we 
can proceed and compute the absolute transverse spectra 
in pA collisions. We assume the proton and the deuteron 
to interact as pointlike objects at an impact parameter b 
with the nucleus. Its nucleons, N, have isospin averaged 
parton distribution functions /^/jv — Z fi/p + {A— Z)fifni 
with A and Z the atomic mass and atomic number. Fur- 
thermore, we assume that j4-nucleus partons scatter only 
once on the proton or the deuteron, due to their small 
density. Then we may generalize Eq. |0J) as follows, with- 
out introducing further free parameters: 



da 



pA^iX 



(PprdycPb 



^ i {^fi/p)yi,PT 



da 



lA 



d'^prdyid'^b 

+TA{b)'^{xfi/A)y,,PT 



(7) 



yi=v 
da'P 



d'^prdyi 



Hadron production is then computed analogously to 
Eq. ®. 

Note that, due to Eq. at large px, or as A-^1 [as- 
suming Ta (6)^5(6)], the ^-integrated pA cross section 
reproduces exactly the pp cross section discussed in the 
previous section. In this way, we can calculate consis- 
tently both the pp and pA transverse spectra in the same 
formalism. 

The Cronin ratio, Rba, of the inclusive differential 
cross sections for proton scattering on two different tar- 
gets, normalized to the respective atomic numbers A and 
B is given by 



RBAiqr) 



B dapA/d'^qrdy 
A daps/d^grdy 



(8) 



First, we can test the GE formalism against low-energy 
data at y/s = 27 A GeV |3| for the ratio of of midrapidity 
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FIG. 2: Cronin effect in charged pion production at = 21 A 
GeV. Plotted is the ratio of the minimum bias charged pions qq^ 
spectrum at r; = in and pBe colHsions. The solid line is for 
a scale choice Qp = = raq-jl and intrinsic (fcy) = 0.52 GeV'^. 
The theoretical error due to the uncertainty in po = O.SitO.l GeV is 
shown as a dotted band. The dashed line shows the result without 
intrinsic kq^ (the theoretical uncertainty is not shown in this case). 
Data points taken from Ref. 0. 



77 = pion spectra in proton-tungsten (jiW) and proton- 
berillium (pBe) collsions: 



R 



■pW/pBe 



(6 = bw) 



dydq^d'^b 



^^pBe~.^^X 

(b = bse), 



dydq^d^b 



(9) 

In our computations we approximated the minimum bias 
cross sections in Eq. by computing Eq. |7J at an 
average impact parameter bw = 5.4 fm and bse = 2.3 
fm, respectively. These values were computed with the 
Monte Carlo model 0|, in order for a pA collision at 
fixed impact parameter to produce the same number 
of participant nucleons as a minimum bias one. The 
result is shown in Fig. [51 The two choices of scale, 
Qp = Qh ~ "niT /2 and Qp = Qh ^ "mr /2 - along with the 
respective fits of po ^^nd K from Table - give approxi- 
mately the same result, and only the former choice is used 
in the figure. The dotted area represents the theoretical 
error due to the uncertainty in the fit of po = 0.7 ± 0.1 
GeV. The computation reproduces satisfactorily the ex- 
perimental data inside the theoretical errors. 

Turning to 7r° production at \/s = 200 GeV at \y\ < 
0.3, in the left panel of Fig.|21wc compare our computa- 
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FIG. 3: Cronin ratio in minimum bias d+Au collisions at y/s = 200 GeV. Left: Cronin effect on neutral pion production. The solid 
line is for Qp = Qf^ = m,y/2 and (kx) = 0.52 GeV^. The theoretical error due to the uncertainty in pQ = 1.0 ± 0.1 GeV is shown as a 
dotted band. The dashed line is computed with no fcy smearing, (k^) = GeV'^, and its theretical uncertainty is not shown. Data points 
are from the PHENIX collaboration, Ref. Illl . Error bars represent statistical errors. The empty bands show systematic errors which 
can vary with qt- The bar at the left indicates the systematic uncertainty in the absolute normalization of the pA cross section. Right: 
Cronin Effect on gluon (dashed line), quark (dot-dashed) and averaged quark and gluon production (solid), with Qp = Qh = mT/2 and 
{k^) = 0.52 GeV2. 



tion for 

^dAu ^ ■ n , ,n, (0 ^OAu) / TaiAOAu) , o , , (10) 

with bAu — 5.7 fm, to experimental data from the 
PHENIX collaboration The results obtained with 

the two choices of scales are similar, and only the com- 
putation with Qp = Qh = rriT/'i- is shown. At this en- 
ergy the sensitivity of the result to the error in the fit 
of Pq and to the scale choice is smaller than at Fermilab 
energy, thanks to the reduced steepness of the pp spec- 
trum. The result is compatible with data on the whole 
Pt range inside the experimental statistical and system- 
atic errors. Despite this caveats, the GE model tends 
to slightly overestimate the data at ^ 2 GeV. What 
we see in Fig. 13 is therefore a possible indication for a 
dynamical shadowing in addition to the basic Glauber 
geometrical shadowing. Its magnitude is consistent with 
the range of dynamical shadowing explored in H, S 
using a variety of shadowing functions [4^ . 

In the right panel of Fig. we plotted the correspond- 
ing Cronin effect at the parton level. The Cronin ratio 
peaks at fairly large ~ 6 GeV, compatible with the 
expected {z) ~ 0.6. The peak in our computation is 
positioned at significantly larger transverse momentum 



than found in the MV model of Ref. Ref. 26]. This dif- 
ference may be due to their choice of parameters: they 
compute the Cronin ratio Rba for a nucleus A such that 
Qs/A = 2 — 3 GeV, and an arbitrary reference B such 
that Qs/B = 1 GeV. Also an infrared cutoff A = 200 
MeV was employed. The same value of the peak [pT — 3 
GeV) was found in the MV model via numerical com- 
putation of Ref. li^l, with a slightly different choice of 
parameters. None of the values of Qs/Ai Qs/b ^^id A 
where fixed by fitting absolute inclusive spectra in pp or 
pA collision. Therefore, the results of Ref. [i^ should 
be understood as illustrative of the qualitative features 
of the MV model. As Fig. O clearly demonstrates, frag- 
mentation strongly distorts the features of the parton 
level Cronin effect. Therefore, the transverse momentum 
scales illustrated in the qualitative saturation models as 
in Refs. [l^ ItH . llsL 126*1 . which up to now do not attempt 
to include the distortions of scales due to hadronization 
processes, should not yet be taken literally. 

To understand better the possible emergence of dy- 
namical shadowing at RHIC, we have two powerful han- 
dles. The first one is the rapidity dependence of the 
Cronin effect, which we will address in a separate publica- 
tion (see also 0, for a discussion in the framework 
of approximated GE, MV and saturation models). The 
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FIG. 4: Central to peripheral ratio for neutral pion production 
in minimum bias d+Au collisions at y/s = 200 GeV. The central 
0-20% bin corresponds to an average b = 3.5 fm and the peripheral 
60-88% bin to an average b = 6.5 fm. The solid line is computed 
with a scale choice Qp = Qi^ = mT/2 The theoretical error due to 
the uncertainty in po = 1.0 it 0.1 GeV is shown as a dotted band. 



second handle is the centraHty dependence of the Cronin 
effect: the more central the collision, the higher the den- 
sity of target partons, the higher the shadowing effects 
induced by non-linear parton interactions. A very nice 
observable will be, in this repect, the ratio of px spec- 
tra in central and peripheral collisions. This ratio has 
the additional experimental advantage that most of the 
systematic errors shown in Fig. |31 are expected to can- 
cel out. This will provide a rather precise comparison to 
the GE model prediction which isolates geometric shad- 
owing only. Our predictions is shown in Fig. ^ for two 
different choices of scale. The average impact parameter 
is b = 3.5 fm and b = 6.5 fm for central and periph- 
eral collisions, corresponding to centrality classes 0-20% 
and 60-88%, respectively [43. In the figure, the theo- 
retical uncertainty due to the uncertainty in the fit of 
Po = 1.0 ± 0.1 GeV is shown as a dotted band. If dy- 
namical shadowing is present, we would expect a larger 
deviation of the plotted curve from the data than what 
is observed in minimum bias collisions in Fig. 13 



CONCLUSIONS 

We have studied the Cronin effect in pA and dA col- 
lisions in the context of Glauber-Eikonal models. These 
models incorporate parton multiple scatterings and uni- 
tarity in pQCD in a consistent way. Moreover, they in- 
clude a detailed parton kinematics and reproduce, in the 
limit of A = 1, the hadron transverse spectra as com- 
puted in the pQCD parton model. The analysis of pp 
spectra allows to fix the free parameter of the model, and 
to compute the spectra in pA collisions without further 
assumptions. 

A powerful feature of GE models is that they automat- 
ically include in the computation geometrical shadowing 
effects induced by unitarity and parton multiple scatter- 
ings. By isolating geometrical shadowing, one can use the 
GE model computations as a baseline in the search for 
genuine dynamical shadowing effects due to non-linear 
parton interactions. 

We tested our computation of the Cronin effect in min- 
imum bias collisions at -^s = 27.4 GeV. The same for- 
malism applied to the recently measured d + Au data 
at = 200 AGeV HI El 111 El describes weU the 
Cronin effect at large pt, with a tendence to overesti- 
mate by '--^ 10 — 20% the effect for tt" at 77 = and 
Pt ^ 2 GeV. Our results are surprisingly similar to pre- 
dictions based on phenomenological approximated GE 
models iSt & li] in spite of the inclusion of geometrical 
shadowing in our GE approach. This provides further ev- 
idence for the possible existence of moderate shadowing 
in the x ~ 0.01 — 0.1 range as explored in those refer- 
ences. However, radical gluon shadowing as predicted 
in is not supported by the data. It remains to be 
seen if the most recent variations of saturation models 
can be fine tuned to account to the thusfar featureless 
RdAu ~ 1 RHIC data. Future analysis of the centrality 
and pseudo-rapidity dependence of the Cronin effect at 
RHIC will provide a powerful tool to further constrain 
the magnitude of the dynamical shadowing effect. 
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